Possible resolution of the Casimir force finite temperature correction "controversies" 
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By considering the effect of diffusion on the external electric field response of charge carriers in 
metals and semiconductors, it is shown that the finite temperature correction proposed Bostrom and 
Sernelius requires substantial modification, and there is no large correction as suggested for good 
conductors. The apparent violation of the Third Law of Thermodynamics of the various proposed 
temperature corrections to the Casimir force is also resolved. Finally, the effect of Debye screening 
on electrostatic calibrations between pure germanium surfaces is calculated. 
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I. INTRODUCTION 



In our recent experimental work aimed at a precision 
measurement of the Casimir force between pure germa- 
nium (Ge) plates Q, we have discovered that the elec- 
trostatic calibration for Ge, a semiconductor, is substan- 
tially different from what was expected, assuming Ge to 
be a good conductor. The problem is that a static electric 
field can propagate a finite distance into a semiconductor; 
this distance is determined by the combined considera- 
tion of diffusion and field driven electric currents, lead- 
ing to an effective field penetration length (Debye-Hiickel 
length) 
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where ct = Ch + Ce is the total carrier concentration, 
which for an intrinsic semiconductor, Ce — Ch. For in- 
trinsic Ge Awl /xm, while for a good conductor, it is less 
than 1 nm. A is independent of the applied field so long 
as the applied field E times A is less than the thermal 
energy, kbT where kb is Boltzmann's constant. In this 
limit, and at sufficiently low frequencies and wavenum- 
bers, thermal diffusion dominates the field penetration 
into the material. A sufficiently low frequency for Ge 
would be Uc/A ~ 10 GHz, where Vc is a typical thermal 
velocity of a carrier. 

An analysis of the electrostatic energy between parallel 
plates is given in the Appendix below, as well as the 
effect of field amplitude on A. This analysis is crucial 
toward our ongoing experimental efforts, especially for 
the electrostatic calibration. 

In light of this analysis, it has become clear that none 
of the recent papers describing finite temperature effects 
on the Casimir force have taken into account the thermal 
diffusion of carriers (electrons and/or holes) in the treat- 
ment of the boundary value problem. As such, a compre- 
hensive review of recent work will not be presented here; 
only the work by Bostrom and Sernelius, that led to the 
recent "controversy," will be discussed 0- 



II. CALCULATION OF THE THERMAL 
CORRECTION 

To calculate the effects of finite temperature, the elec- 
tromagnetic mode photon excitation number of 1/2 due 
to zero point fluctuations is replaced by 



n{ijj) = coth 
which has simple poles at 
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where in the following discussion we take only integers 
n > 0. 

Following [^], the integral over oj in determining the 
field energy between two flat plates is replaced by the sum 
over the poles that occur at the Matsubara frequencies. 



£ = kbTj2 /K) 
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where the prime indicates a factor of 1/2 for the n = 
term, and for the TE modes (electric field parallel to the 
surface) 



(5) 



where q represents the electromagnetic field wavenumber 
in the space between the plates, in direction perpendic- 
ular to the plates. 

The function in the integral, for the TE mode, is 



G^^(c.„,,) = l-(:^] e-^-o. 
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where i = 0, 1 represents the space between the plates (0) 
or inside the plates, (1), and ti is the respective electric 
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permittivity along the imaginary frequency axis, and d is 
the plate separation. 

It is argued in [3| that for realistic materials, 
G^^{0,q) — 1 and hence does not contribute to the en- 
ergy that leads to the Casimir force between the plates. 
For example, the low-frequency permittivity of a metal 
is given by the conductivity cr. 
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for which 71 = 72 at = O.For distances greater than 
a few /i m, for T = 300 K, the net force is reduced by 
a factor of two compared to what is expected for a near 
perfect conductor if both the TE and TM modes are 
included. 

This result is in contradiction to experimental results, 
particularly There has been much discussion of this 
correction, but to now, the effect of diffusive field screen- 
ing has not be taken into account in a satisfactory man- 
ner, or at all. 



III. INCLUSION OF THE DEBYE SCREENING 
LENGTH 

For realistic conducting materials (metal, semiconduc- 
tor), Eq. ([7|) is not possibly correct, for we know an 
electric field near the surface causes charges to move, 
and this tends to screen out the applied field. Electric 
fields applied either parallel or perpendicular to the sur- 
face will be screened, varying compared to the field Eq 
at the surface, in the material, as 



E{x) = Eae 
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which is valid wheni?oA < k^T (or EqX < k^T/ei for 
fields perpendicular to the surface). 

The fact that electric fields do not penetrate any appre- 
ciable distance into even very poor conductors is exper- 
imentally well-known. Eq. ([7|) incorrectly describes the 
observed penetration of low frequency fields into conduc- 
tors. 

The screening effect can be treated in a heuristic fash- 
ion. The exact solution to the combined electrodynamic 
and diffusion problem is beyond the scope of the present 
note, the intent of which is to illustrate an effect that 
has been overlooked. Noting that A is independent of 
frequency at low frequencies (< 10^° Hz), and that E 
and its derivative are continuous across the boundary, 
we can rewrite Eq. [7]) to include the spatial variation 
due to Debye screening. 



(10) 



assures continuity of the derivative across the boundary. 
The point is that the wavenumbcr in the material cannot 
be less than 1/A. 

For this modified function, for ^ 0, G"^^ is equiva- 
lent to the perfect conductor result. Thus the large cor- 
rection suggested in |3] applies only to materials where 



charges are not free to move, and diffusive effects do not 
enter. We can question whether this condition can ever 
really be met, but for any realistic slightly conducting 
material at finite temperature, there will always be a fi- 
nite screening length, and hence a full contribution from 
the G^'^, n ~ Q mode, for good conductors where A is 
large to frequencies of order 10^^ Hz. A detailed calcu- 
lation for Ge is required because A starts falling off for 
frequencies above 10 GHz, in the region where the ther- 
mally excited photons contribute most significantly 



IV. VIOLATION OF THE THIRD LAW? 

It has been suggested that the w = term in Eq. () 
show a manifest violation of the Third Law of Thermo- 
dynamics (sometimes referred to as the Nernst Heat The- 
orem) because the system entropy (identifying £ as the 
free energy F), given by 
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is not zero in the limit T —t O.This conundrum has been 
addressed and clarified in 5], but there is a simpler ar- 
gument that will be presented here. 

In particular, we can question what it means to convert 
the integral over uj into a contour integral and hence sum 
over the Matsurba frequencies. If we did not convert this 
integral into a sum, the separate identification of the zero 
frequency contribution would not be made. Specifically, 
T docs not appear in the Casimir force calculation expect 
through the mode photon number. Differentiating Eq. 
(2) with r, we find 
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which indeed goes to zero at T = 0. 

Alternatively, if we consider the entire sum, Eq. (4), 
let T go to zero, we find. 



/cfcT^ /(w„) = fcfcT / f{2TTnkbT/h)dn 

n=0 
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where the simple substitution 27rnfcf,T/fi, = lj was made. 
Hence, the temperature does not appear explicitly in the 
total free energy, and the entropy indeed goes to zero at 
zero temperature. 

The apparent violation of the Third Law is due simply 
to the isolation of a single term in the total free energy. 
Not considering the entire system in calculating the en- 
tropy is generally considered a sophomoric error. 



V. CONCLUSION 

By including the effect of charge movement and screen- 
ing through the Debye length, it is shown that the large 
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correction to the Casimir force predicted in is not ap- 
plicable to realistic materials. It should be noted that 
these correction apply to all conductors when the dis- 
tance scale approaches the Debye length, which for a 
good conductor is 0.1 nm. The Debye length is constant 
in good conductors up to frequencies of order 10^"' Hz, so 
we can expect the full perfectly conducting force for any 
metal, at large separations. 

The case of a semiconductor like Ge is slightly more 
complicated because A begins to increase for frequen- 
cies of order 10 GHz, so a detailed anlysis of this case 
is required. However, it can be expected that there is a 
significant contribution from the TE, n = mode. 

It is also shown that the apparent violation of the 3rd 
Law of Thermodynamics is due to the isolation of a single 
term in an expansion that becomes an integral in the 
limit of T ^ 0. 



VI. APPENDIX 

The potential in a plane semiconductor, if the potential 
is defined on a surface x = is 



V{x) = F(0)e-l^l/^ 



(14) 



where A is the Debye-Hiickle screening length, defined 
previously. 

We are interested in finding the energy between two 
thick Ge plates separated by a distance d, with a voltages 
+V/2 and —V/2 applied to the backs of the plates. In 
this case, the field is normal to the surface. After we find 
the energy per unit area, we can use the Proximity Force 
Theorem to get the attractive force between a spherical 
and flat plate. 

Let a; = refer to the surface of the plate 1, and x — d 
refer to the surface of plate 2. By symmetry, the potential 
at the center position between the plates is zero. The 
potential in plate 1 can be written 



^1(2:) = V/2 - {V/2 - V;)e-I^l/^ 
and the space between the plates 



(15) 



^0(2;) = -2Vsx/d+Vs 

where we assume the field is uniform. Vs, the surface 
potential, is to be determined. 



We need only consider the boundary conditions in 
plate 1, which are 
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1/0(0) = Fi(0) 
(which has already been used) 

dVi{x). _dVo{x)^ 
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ax, ax 

where the last two imply that D 
across the boundary. 
The solution is 



\x=0 

eE is continuous 
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With this result, it is straightforward to calculate the 
total field energy per unit area in both plates and in the 
space between the plates. The result is 
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where the dimensionless length y — ed/X has been intro- 
duced. 

If V — Vg is large compared to fcf,T, the effective pen- 
etration depth increases because the charge density is 
modified in the vicinity of the surface. The potential in 
the plates is no longer a simple exponential, however one 
can define an effective shielding length [6|] 
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where 



V/2 - Vs 
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with the results plotted in Fig. 2. Given that k^T = 30 
meV, at plate separations of order 1 /xm for Ge this begins 
to be a large correction when voltages larger than 60 mV 
are applied between the plates at separations of order 
1 /xm. 
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FIG. 1: The separation determined by the force resulting from 
an applied potential, assuming a 1/d force variation, should 
be multiplied by this value. For Ge, ed/A = 20 corresponds 
to 1 fim plate separation. 



FIG. 2: Effective field propagation distance as a function of 
the dimensionless potential difference between the bulk and 
surface. 



